Neutrino orbital angular momentum in a plasma vortex 
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It is shown that an electron-neutrino beam, propagating in a background plasma, can be decomposed into 
orbital momentum (0AM) states, similar to the OAM photon states. Coupling between different 0AM neutrino 
states, in the presence of a plasma vortex, is considered. We show that plasma vorticity can be transfered to the 
neutrino beam, which is relevant to the understanding of the neutrino sources in astrophysics. Observation of 
neutrino OAM states could eventually become possible. 

PACS numbers: 



I. INTRODUCTION 



Neutrino beam interaction with a dense plasma has received 
considerable attention in recent years. In particular, it was rec- 
ognized that collective neutrino-plasma displays strong simi- 
larities with laser-plasma interaction phenomena. This is due 
to the fact that neutrino dispersion relation in a plasma is for- 
mally analogous to the photon dispersion relation [1]. Cou- 
pling between neutrino beams and dense plasmas can then be 
conceived HUH, which is particularly relevant for the under- 
standing of supernova type II explosions, where an intense 
and short neutrino burst is emitted [4]. It was also recognized 
that neutrino ponderomotive force can excite kinetic plasma 
instabilities Jj, @], resulting from negative neutrino Landau 
damping of electron plasma waves, similar to that observed 
for photons [7]. Furthermore, it is also possible to show that 
neutrinos acquire an effective electric charge in a plasma, in 
a way very similar to photons JH @]. Another important as- 
pect is related to the possible excitation of relativistic plasma 
wake fields by a neutrino burst ITolfTlll . and the generation of 
magnetic fields [12], through processes that mimic those of 
photon beams in a plasma. 

In this work, the similarities between neutrino and photon 
dispersion are explored even further, by introducing the con- 
cept of neutrino orbital angular momentum (OAM). This is a 
new aspect of neutrino plasma physics, which is related with 
the possible exchange of orbital angular momentum between a 
neutrino beam and a rotating plasma. It is known that photons 
can carry, not only intrinsic angular momentum or spin, which 
is associated with their polarization state, but also orbital an- 
gular momentum (OAM) 11311 . The existence of OAM pho- 
ton states has been experimentally demonstrated early in 1936 
Jl4ll . but only recently photon OAM momentum received at- 
tention, after the demonstration that Laguerre-Guassian laser 
modes correspond to well defined OAM modes, and that these 
photon modes not only can be measured as a photon beam 
property lfl5l [l6ll . but can also be detected at the single pho- 



ton level [ 1" 

Utilization of photon OAM in the low frequency radio wave 
domain was recently proposed by one of us [18], as an ad- 
ditional method for studying the properties of radio sources. 
Here we show that a neutrino beam can be decomposed into 
OAM momentum states, similar to the OAM photon states, 
and we study the coupling between these modes, when the 
neutrino beam propagates across a plasma vortex region. We 
will also show that plasma vorticity can be transferred to the 
neutrino beam, which can likewise be of great relevance to the 
understanding of the astrophysical sources of neutrino emis- 
sion. Detection of finite neutrino OAM states would give us 
additional information on the collapsing star, in the case of 
a supernova explosion, and would allow us to measure the 
vortex properties of the surrounding dense plasma. Observa- 
tion of neutrino OAM states could eventually be possible, by 
studying the conservation of total angular momentum at de- 
tection, and by adapting the existing detection schemes to that 
purpose. Detection of radio signals from neutrinos interacting 
with the Moon, as recently proposed lfl9ll . could provide one 
possible method for measur ing such neutrino OAM states. 



n. BASIC EQUATIONS 

We know that the spinor field y/ describing electron- 
neutrinos moving in a dense plasma can be described by the 
Dirac equation 



ih d~t W = HV 



H = j5m v c 2 + a • pc + V 



(1) 
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where m v is a non-zero rest mass, j3 and a are the well known 
operators. For simplicity, the plasma medium is assumed 
isotropic, and magnetic field effects are ignored. We retain 
explicit physical units, with c 7^ 1 and hj^l, for future com- 
parison between neutrino and photon states in a plasma. The 
effective plasma potential V(r,f), results from the weak inter- 
actions between the neutrinos and the background electrons 
and ions (or protons) of the background plasma. We know 
that, for an electrically neutral medium, the proton contri- 
butions are exactly canceled by the neutral part of the elec- 
tron contribution, and the resulting plasma potential is simply 
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*iven by 

V = V2G F n e (r,t)=hg(0 2 p 



= V2 



G F (2) 



where n e is the electron mean plasma density, co p is the 
electron plasma frequency, e and m e are the electron charge 
and mass, Eq is the permittivity of vacuum, and Gf is the 
Fermi constant. Equations (Q]l-© have been used to study the 
collective decay of plasmons into neutrino-antineutrino pairs 
Il20l Bill . Here, we neglect such neutrino-antineutrino cou- 
pling processes. We also neglect mass oscillations, but their 
possible contribution will be discussed later. For a constant 
plasma density, we can use the Foldy-Wouthuysen procedure, 
leading to a diagonalized Hamiltonian 



H = j3 



- p 2 c 2 - 



■V 



(3) 



We will focus on a single spin (or helicity, in the limit of a 
negligible mass) electron-neutrino state. The corresponding 
solution can then be written as yr(r,t) — I'oexp^'k ■ r — i(Ot), 
where the neutrino frequency co=<H > /ft, and the wave vec- 
tor k =< p > I ft satisfy the neutrino dispersion relation 



(CO 



-gco 2 p ) 2 = k 2 c 2 - 



col 



(4) 



with <»o = m v c 2 /ft. This shows a strong similarity with the 
photon dispersion relation in a plasma. Let us now consider 
a non-homogeneous and non-stationary medium, where the 
plasma density and the plasma potential are no longer con- 
stant, but vary on time a length scales much longer than the 
neutrino period and wavelength. 



V(r,f)=Vo + V(r,f) 



(5) 



which is the consequence of a space-time variation of the 
background plasma density and of the corresponding plasma 
frequency, as described by 



fflp = «#)[l+e(r,f)] = fl> 



1 , X 

H R(r,t) 



(6) 



where n is the electron density perturbation. By comparing 
the two expressions, we conclude that 



V(r,t) = ftgco 2 p0 e(r,t) = hg^n(r,t) (7) 

n 

In order to solve the neutrino field equation, we can use the 
quasi-classical or WKB approximation and make the replace- 
ments 



H -> ft { co- 



'dt 



ft(k+iV). 



(8) 



co- /-J = /3^a) 2 + (k + /V) 2 c 2 + -(V + V) ( 9 ) 



The Hamiltonian operator (f3]l then becomes 



Assuming that the dispersion relation is always locally sat- 
isfied for the non-perturbed plasma potential Vo, we can use a 
WKB solution of the form 



V(r,f) = v F(r,f)exp(/k-r-/o)f) 



(10) 



where *F(r,?) is now a slowly varying amplitude. Expanding 
the operator (O we can then write the evolution equation for 
this amplitude as 



d 

ift—^V 
dt 



(co-gcof, 



Zk-V-^V 2 



c 2 (k-V) 2 
2{co-gcol) 



■V 



(11) 

Noting that the neutrino mean velocity can be identified with 
the group velocity associated with the dispersion relation ©, 
we get 



c 2 k 



v v = 



dco 



<9k co-gcoj 



(12) 



We can also introduce the neutrino relativistic factor y v , as 
defined by 



m v y v c 2 = ft(co-gco 2 ) 



(13) 



Using these quantities, we can rewrite equation dTTb in a more 
appropriate form 



ft 



2m v y v 



v 2 + (^-v) 2 + ^„ 

\ c ) no 



(14) 

In the limit of a negligible neutrino velocity, this would take 
the form of a Schroedinger equation. A similar equation was 
also derived for the case of a quantum free electron laser I22I 
I23IL Let us first assume that a neutrino beam propagates along 
a given axis Oz, in the absence of any plasma perturbation, n = 
0. Retaining only the first order derivatives in the z- direction, 
we can then reduce the wave equation to a paraxial equation, 
similar to that describing a laser beam in the focal region 



2ik^ ) m, = o 



(15) 



where the wavenumber k obeys the dispersion relation ©. We 
can try a Gaussian solution of the form 



y/b(r) = VoO,z) = *Fo(z)exp 



-2 1 



ik 



2R(z) 



(16) 



Using the transverse Laplacian in cylindrical coordinates, we 
realize that this solution satisfies the paraxial equation, if we 
assume that 



R(z)=Ro + (z-zo) , *Po(z)=mo 



(17) 



where Rq and uq are constants. Very simple arguments can be 
used to justify the relevance of such Gaussian beam solutions 
in neutrino physics. For instance, the intense neutrino burst 
emitted by a collapsing type II supernova, is naturally focused 
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by the plasma inhomogeneities of the star material. The neu- 
trino dispersion relation (|4]i shows that neutrinos will be fo- 
cused into the lower density plasma regions, in the same way 
as high energy photons would do. Another reason for neutrino 
focusing into a nearly Gaussian form, would be the occurrence 
of filamentational instability of the ultra intense neutrino burst 
coming out of the collapsing star. A more general cylindrical 
type of solution of the paraxial equation ( TT3b can be repre- 
sented in the basis of orthogonal Laguerre-Gaussian modes, 
as defined by 



V(>V)=E <t>piF p i(r,(j>)exp(ikz-i(Dt) 
P i 



(18) 



with 



'>'"••* ' X (5)"^'(5) eXp( ^"2^ ) (19) 



The relevance of this type of representation will become ap- 
parent below. We recover the particular Gaussian solution by 
taking / = 0, and R — —iw 2 . The new quantity w defines the 
neutrino beam waist. By making an appropriate choice of the 
normalization factors multiplying equation ( fT9l . we can verify 
the orthogonality conditions 



> r2n 

rdr d<$>F* pl (r, <j) )F p>v (r, 0) = 8 pp , 8„, 



(20) 



III. NEUTRINO MODE COUPLING 

Let us now consider the case where the electron density and 
the plasma potential are not constant, n^O. In cylindrical 
geometry, we can use a general representation in terms of the 
Laguerre-Gaussian modes, as 

n(r,*)=£ Jj^J ^ni(r,q,Q.)exp{il(t>+iqz-i£lt) (21) 

To be more specific, we can concentrate on helical type of 
perturbations described by 



n(r,f) = n(r)exp ^— ^2 J cos ('o0 + "? z — '^0 (22) 

In the presence of such density perturbations, the different 
Laguerre-Gaussian modes of the neutrino field will be cou- 
pled, and the generic wavefunction solution will take the form 



V(r,f) = £ yr np i(r)exp[iS(r,t) 

n,p,l 



(23) 



where y/„ p / are slowly varying ampitudes, and S(r,t) is the 
phase function as defined by S(r,t) = k„z—(O n t. Here we have 
used co„ = CO + n£l, and the corresponding wavenumbers k„ 
are determined by the neutrino dispersion, with CO replaced by 
CO,,, valid for the unperturbed plasma density no. We can use, 
for the above amplitudes 



Wn P i(r) =^n P i(z) F p i(r,ij>) 



(24) 



Replacing this in the paraxial equation, where the term corre- 
sponding to the perturbed potential is retained, we obtain the 
following mode coupling equations 

d 



n'p'l' 



(25) 



where we have defined the wavenumber mismatch A nn i 
(k„i —k n ±q), and used the following coupling coefficients 



K(npl,n'p'l') = 



gCOn <*>p0 

k„c 2 n 



8{n-n'±\) 



2 IT 



dcj> I rdrn^F^^F^sxpliil^lo-l)^} (26) 
/o Jo 1 

Here, in the expressions of F* { (r) and F p ni(r) we have taken 

out the factors exp(— il<j>) and exp(il'<p). The coupled mode 

equation ( |25T l shows that, if we start from an initial neutrino 

state with no orbital angular momentum, /' = 0, we will excite 

states of finite angular momentum / ^ 0, due to the existence 

of a perturbed plasma helical structure. But such coupling has 

to verify certain resonant conditions. Depending on the radial 

profile of the plasma vortex «(r), various modes p^p' can be 

excited. But we take here the simplest case of p — p' ', which 

allows us to drop this index from the coupled mode equations, 

and write them in a much simpler form 

d 



dz 



M»v(z) 



iK 



^ v+1 (z)e ,az + W v ^(z)e 



-iAz 



(27) 



Here v is an integer, and we have used a simplified notation 
for the amplitudes, ^Py = m n+v „ i+ v u . Approximate expres- 
sions for the coupling coefficients and wavenumber mismatch 
can be written as 



K 



gco 



CO 



kc z no 



lp0 »(o) 



Q 



(28) 



For an exact phase matching, such that q = —£l/c, the cou- 
pled mode equations d2Tb can easily be integrated in terms 
of Bessel functions. For initial conditions such that ^/(O) = 
^Po^/o, corresponding to an initial Gaussian neutrino beam, we 
get the following amplitude for the various interacting modes 



>F v (z)=/ Vl IVv(2 J ft) 



(29) 



This shows that the states of higher neutrino OAM can be pop- 
ulated by a linear cascading process, approximately described 
by this simple law. Such a simple but physically relevant solu- 
tion is valid for arbitrary interacting distances in the case of a 
perfect phase matching A = 0. And it remains valid for a finite 
phase mismatch, for interaction distances much shorter that 
2n/A. However, in this more general case and for arbitrarily 
long distances, the transfer of OAM from the plasma rotat- 
ing vortex to the neutrino beam will eventually be blocked by 
phase mixing. 



IV. CONCLUSIONS 

In this work we have considered the orbital angular mo- 
mentum of a neutrino beam in a plasma, to our knowledge for 
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the first time. We have studied the interaction of the neutrino 
beam with a plasma vortex, and explored the similarities be- 
tween neutrino and photon dispersion in a plasma. We have 
shown that a neutrino beam, propagating in a dense plasma 
background, can be generally described by a superposition of 
OAM states, similar to the photon OAM states recently ex- 
plored in the literature. The existence of plasma vortices in- 
troduces coupling between different neutrino OAM states. We 
have shown that plasma vorticity can be transfered to the neu- 
trino beam. This process can be of particular relevance to the 
understanding of the astrophysical neutrino sources. Study 
of total angular momentum conservation, using the existing 
neutrino detection systems, could eventually lead to the first 



observation of neutrino OAM states. 

In this work, a simplified neutrino description was used. 
In particular, we have neglected the neutrino mass oscillation 
process. It is known that in a dense plasma, these oscilla- 
tions are resonantly enhanced by the so called MSW effect. 
This means that the existence of plasma vorticity in the res- 
onant plasma region where the MSW effect takes place will 
be able to excite OAM states of the different neutrino species. 
Generalization of the MSW effect, taking the OAM states into 
account, is in progress and will be the subject of a future work. 
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